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Abstract. We investigate the relation between the existence and size of radio halos, which are believed to be 
created by star formation (SF) related energy input into the interstellar medium, and other galaxy properties, 
most importantly star formation activity and galaxy mass. Based on radio continuum and Ha observations of a 
sample of seven late- type spiral galaxies we find a direct, linear correlation of the radial extent of gaseous halos on 
the size of the actively star-forming parts of the galaxy disks. Data of a larger sample of 22 galaxies indicate that 
the threshold energy input rate into the disk ISM per unit surface area for the creation of a gaseous halo depends 
on the mass surface density of the galaxy, in the sense that a higher threshold must be surpassed for galaxies with 
a higher surface density. Because of the good prediction of the existence of a radio halo from these two parameters, 
we conclude that they are important, albeit not the only contributors. The compactness of the SF-related energy 
input is also found to be a relevant factor. Galaxies with relatively compact SF distributions are more likely to 
have gaseous halos than others with more widespread SF activity. These results quantify the so-called "break-out" 
condition for matter to escape from galaxy disks, as used in all current models of the interstellar medium and 
first defined by Norman & Ikeuchi (1989). 
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> 1. Introduction 
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The investigation of gaseous halos in late-type spiral 
galaxies is closely linked to studies of the galaxies' chemi- 
' cal evolution. Halos are potential mediators of metal redis- 
tribution via disk- halo interactions (e.g. Bregman 1980). 
If escape velocity is reached, for example in galactic winds, 
outflowing halo gas influences the evolution of the inter- 
galactic medium by means of metal and energy injection, 
thus adding to the potential importance of gaseous halos 
on cosmological scales (e.g. Heckman 2005). 

Over the last decade it has become clear that spiral 
galaxies with high-mass star formation (SF) have gaseous 
halos (e.g., Heckman, Armus, Miley 1990; Lehnert and 
Heckman 1995; Dahlem 1997; Rossa and Dettmar 2003a). 
In fact, we have found evidence (Dahlem et al. 1998; 
Dahlem et al. 2001) that almost all far-infrared (FIR) 
"warm" edge-on spirals (i.e., with 60 /in to 100 (im FIR 
flux ratios feo/fwo ~ 0.4) have gaseous halos. The only 
remaining exceptions from this rule are galaxies for which 
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no firm statement can be made due to the lack of data of 
sufficient quality. 

These gaseous halos comprise all known phases of the 
interstellar medium (ISM) previously found in the disks 
(e.g., Dahlem 1997). Optical (Ha) emission was detected 
in M82 already by Lynds and Sandage (1963). Radio de- 
tections came later, starting with NGC4631 (Ekers and 
Sancisi 1977). Yet later, X-ray emission was detected from 
halos (e.g., Watson et al. 1984). 

Further progress in this field of research has been 
made by studies in various wavebands, such as optical 
Ha imagery (Rossa and Dettmar 2000, 2003a, b; see also 
the work by Dettmar 1992 and Rand 1996) and X-ray 
imaging spectroscopy (e.g., Dahlem et al. 1998, Pietsch et 
al. 2000, Strickland et al. 2004a, b and references therein; 
Tullmann et al. 2006). HI was detected in halos around 
spirals (e.g., Fraternali et al. 2004 and references therein; 
Boomsma et al. 2005) and, most recently, diffuse ultra- 
violet (UV) emission was found from the starburst out- 
flows of NGC253 and M82 (Hoopes et al. 2005). 

Radio synchrotron emission from relativistic electrons 
is the most extended and pervasive component of these ha- 
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Table 1. Galaxy Sample 





NGC 1808 


M82 


NGC 4666 


NGC 4700 


NGC 7090 


a (2000) 


05:07:42.3 


09:55:52.2 


12:45:08.3 


12:49:07.6 


21:36:28.6 


5(2000) 


-37:30:47 


+69:40:47 


-00:27:51 


-11:24:47 


-54:33:24 


t>hci (km s" 1 ) 


1000 


203 


1520 


1404 


849 


D (Mpc) 


10.9 


3.2 


26.4 


25.5 


11.7 


PA (°) 


320 


60 


225 


48 


130 




60 


90 


78 


90 


90 



The entries in the rows of Tableware: 

q, <5(2000): J2000 equinox centre coordinates 

Mhei: Heliocentric velocity, from Koribalski et al. (1993; NGC 1808), de Vaucouleurs et al. (1991; M82), Walter et al. (2004; 
NGC 4666) and Dahlem et al. (2005; NGC 4700 and NGC 7090) 

D: Distances are based on H = 75 km s _1 Mpc -1 and a virgocentric infall velocity of 300 km s _1 (cf. Dahlem et al. 2001). 
i: Inclination angle 



los, which makes them most easily detectable. Thus, recent 
attempts to detect radio emission from extraplanar gas in 
actively star-forming galaxies have been quite successful 
(e.g. Irwin et al. 1999; Dahlem et al. 2001). Nevertheless, a 
debate is still continuing about how such halos are created 
and energetically maintained (Dahlem 1997). In a first at- 
tempt to establish general rules of behaviour for gaseous 
halos of spirals, we (Dahlem, Lisenfeld, Golla 1995; here- 
after DLG95) investigated the dependence of the proper- 
ties of radio halos on the level of SF in the underlying 
disks of NGC 891 and NGC 4631. At the time, these were 
the only two galaxies for which we had data of sufficient 
quality to conduct such studies. 1.49 GHz radio continuum 
images were used to measure the properties of the radio 
halos of both galaxies. The most suitable tool to deter- 
mine the radius inside which high-mass SF is occurring, 
rsF, are FIR data. However, only very few images with 
sufficient angular resolution exist (in particular those ob- 
tained with the IRAS CPC detector; see van Driel et al. 
1993). In the case of NGC 4631 we substituted the FIR 
data with an Ha image. 

As part of a small radio continuum survey of FIR- warm 
edge-on spiral galaxies, radio images of several more suit- 
able galaxies were obtained (Dahlem et al. 2001, 2005). 
Here we present an analysis of those radio images with 
sufficient resolution for detailed studies of halo properties. 
We compare these with Ha emission line images and/or 
other tracers of SF in galaxy disks, to investigate whether 
these galaxies show a similar dependence of their radio 
halo properties on the level of activity in the underlying 
disks as NGC 891 and NGC 4631 (DLG95). To ensure that 
processes related to high-mass SF dominate the energy 
balance of the ISM in these galaxies, systems with known 
luminous active galactic nuclei (AGNs) and closely inter- 
acting galaxies were excluded (see Dahlem et al. 2001). 

2. Database 

Here we concentrate on radio continuum and Ha ha- 
los, but different gaseous phases are normally found as- 



sociated with each other (Dahlem 1997; Veilleux et al. 
2005; Tiillmann et al. 2006). Therefore, it is very likely 
that galaxies with either radio continuum or Ha halos in 
general have multi-phase gaseous halos, including all the 
phases of the ISM mentioned above. 

2.1. Sample selection 

Only galaxies with known extraplanar emission both in 
the radio continuum and optical emission lines were con- 
sidered. Also, these galaxies were selected to be neither 
closely interacting (and thus not tidally disturbed) nor to 
host a dominant AGN and are therefore suitable for study- 
ing the dependence of gaseous halo properties on the SF 
activity in the underlying disks. The available datasets 
must allow the measurement of at least the following 
three quantities: 1. The maximum radius at which a spiral 
galaxy actively forms stars, rgp; 2. the maximum radius 
at which nonthermal radio continuum emission from cos- 
mic rays is observed in their disks, tcr; 3. the maximum 
radius at which synchtrotron radio continuum emission is 
detected in the halo, ;*h a i . These selection criteria severely 
limit the number of suitable targets. At present the list of 
galaxies fulfilling all selection criteria, for which we have 
suitable datasets (in addition to NGC 891 and NGC 4631; 
DLG95), comprises five new objects; their basic properties 
are collected in Tabled For these galaxies the following 
observational datasets are used. All data are available to 
us in electronic, reduced form, except the radio continuum 
data for NGC 1808 by Saikia et al. (1990). For this galaxy 
the derived radii were taken from the literature. 

2.2. Radio continuum observations 

Low- frequency (1.4 GHz) radio continuum observations 
are used to trace cosmic-ray (CR) electrons in halos 
around the disks of the sample galaxies. 

The radio data analysed here have all been published 
before. The radio data of NGC 1808 were published by 
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Fig. 1. VLA 1.49 GHz radio continuum image overlaid on 
our Ha image of NGC4700. Contours range from -0.2, - 
0.14, 0.14 (2.5-cr), 0.2, 0.28, 1.6 mJy beam" 1 , spaced by 
factors of \/2, while the logarithmic Ha grey scales range 
from to the maximum observed surface brightness. The 
resolution (full width at half maximum) of the radio image 
is 13" (1.6 kpc). 



Saikia et al. (1990), the data of M82 (=NGC3034) were 
obtained by Reuter et al. (1992). The radio image of 
NGC 4666 is from Dahlem et al. (1997). The 1.4 GHz maps 
of NGC 4700 and NGC 7090 are from Dahlem et al. (2001). 
The observations and their calibration and data reduction 
are all described in the respective papers. 

2.3. Ha images 

Ha imagery is used as a tracer of active SF in the disks 
of galaxies. Almost all Ha images used for comparison 
with the radio continuum data have been taken from pre- 
vious studies. The Ha image of NGC 1808 was published 
by Jimenez-Bailon et al. (2005). M. Lehnert kindly pro- 
vided us with a wide-field Ha image of M 82 (Lehnert et 
al. 1999). The Ha frame of NGC 3175 is from Ryder and 
Dopita (1993). The NGC 4666 Ha image was originally 
published by Lehnert and Heckman (1995) and has al- 
ready been used by us earlier (Dahlem et al. 1997). The Ha 
images of NGC 4700 and NGC 7090 originate from Rossa 
and Dettmar (2003a,b). 

2.4. Additional data 

In addition we include some galaxies in the present study 
for which the measurement of the three radii ri ia io, ^sf 
and rcR (see Sec. I2.1|) is not possible from the available 
images. Still, observations of these galaxies can be useful 
and are added here (once the existence of a gaseous halo 




RIGHT ASCENSION (J2000) 



Fig. 2. Placement of a 12"-wide cut through the Ha emis- 
sion distribution of NGC 4666 along the disk plane, indi- 
cating the area over which the cut displayed in Fig. was 
determined. See Table [21 for details. 



is established) to investigate the dependence of their halo 
properties on some global galaxy parameters, such as e.g. 
their average energy input rates and galaxy mass. 

To this end we use additional 1.4 GHz radio contin- 
uum images of NGC 1406, NGC 3175 and NGC 7462 from 
Dahlem et al. (2001). Furthermore data of three galaxies 
from the sample by Irwin et al. (1999) are used. Their anal- 
ysis of the radio emission is different from ours (Dahlem 
et al. 2001) and not easily comparable. Therefore, we only 
include here galaxies from their sample not showing evi- 
dence of emission beyond the modelled thin disks, which 
is a robust criterion similar to the one used by us in the 
rest of the sample. An Ha image of NGC 7462 from Rossa 
and Dettmar (2003a, b) was also used. 

3. Results 

In the publications from which the data were taken most 
of the basic emission properties of the sample galaxies are 
listed (e.g., total flux densities, centre positions, etc.). An 
example of a (previously unpublished) dataset used here 
is provided in Fig. ^ where we display a contour over- 
lay of a new VLA 1.49 GHz continuum map of NGC 4700 
on our Ha line image. One can see directly related radio 
continuum and optical line emission, not only in the disk 
plane, but also beyond. Other galaxies in our sample ex- 
hibit a similar strong association of radio continuum and 
Ha emission. 

In the following sections we report on our measure- 
ments of the radial extent of both the star-forming parts 
of the galaxy disks, rgF, and of the associated radio syn- 
chrotron halos, rhaio- A correspondence between these two 
radii is a crucial test of the dependence of the existence 



4 



M. Dahlem et al.: Properties of radio halos 



Table 2. Parameters for Cuts 





(I) 


(2) 




(4) 


V J J 






f7) 




Object 


Band 


Location 


HPBW 


Width 


Measurements 




z-offset 








Disk/Halo 


( ") 


( ") 


Averaged 




(HPBW) 


fkDc") 


NGC 1808 


1.4GHz 


Disk 


0.6 


5 


2 













1.4GHz 


Halo 


1 


5 


2 


10 


10 


0.5 




Ha 


Disk 


1 


5 


2 











M 82 


1.4GHz 


Disk 


3.8 


5 


2 













1.4GHz 


Halo 


3.8 


5 


2 


30 


7.89 


0.5 




Ha 


Disk 


1 


5 


2 











NGC 4666 


1.4GHz 


Disk 


14 


12 


2 













1.4GHz 


Halo 


14 


12 


4 


32 


2.29 


4.1 




Ha 


Disk 


1.9 


12 


2 











NGC 4700 


1.4GHz 


Disk 


13 


12 


2 













1.4GHz 


Halo 


13 


12 


3 


16 


1.23 


2.0 




Ha 


Disk 


1 


12 


2 











NGC 7090 


1.4GHz 


Disk 


14 


30 


2 













1.4GHz 


Halo 


14 


30 


3 


60 


1.74 


3.4 




Ha 


Disk 


1 


30 


2 












The entries in the columns of Table |5| are: 

(1) Object name; (2) Observing waveband; (3) Location of cut(s) from which measurements were obtained; (4) Angular resolution 
of data used; (5) Width of the cut over which datapoints were averaged; (6) Number of measurements used to derive mean 
value; (7) Offset of cuts from galaxy disks, in units of arcsec, half power beamwidths and kpc. 



of radio halos on the level of SF in the underlying disks. 
At the same time we measure the radial extent of the ra- 
dio continuum-emitting disks of the galaxies, rcR, and its 
difference compared to the size of the star- forming part, 
Ar = r CR - r SF . 

The measurements are performed by creating slices 
along the disk planes of the sample galaxies to deter- 
mine and and r§F, and additional slices parallel to 
the disk to determine rj la i . The slices have a narrow 
width (approximately one beamwidth of our radio contin- 
uum images), over which surface brightnesses are averaged 
and then plotted as a function of galactocentric distance. 
Details on the placement and size of the regions are sum- 
marised in Table |21 The procedures for measuring these 
quantitites are exactly the same as used by us in DLG95. 
In the following subsections we will use data of NGC 4666 
as an example to demonstrate these techniques. 

3.1. Size of the actively star-forming vs. the radio 
continuum-emitting parts of galaxy disks 

In principle, the maximum radii at which SF is observed, 
rsF, can be measured from our radio continuum maps, at 
the radial points where the surface brightness of the disk 
emission drops drastically. However, in order to make sure 
that these measurements are truely independent of those 
of the radial extent of the radio halos (below), we choose 
instead to use Ha emission from H II regions in the disk as 
a tracer of SF and thus as a measure of tsf- A quantitative 
measurement can be obtained by setting a threshold at a 
given Ha surface brightness, which can be converted into 
a SF rate. This would require the use of flux-calibrated 



Ha images, which we do not have available for all tar- 
gets at this time. The placement of the slice through our 
Ha image along the disk plane of NGC 4666 is displayed 
in Fig. |21 Fig. |3| shows the resulting radial profile. Two 
measurements, one on each side of the disk, are obtained, 
which are then averaged. The cutoff criterion for where the 
star-forming disk ends has been set to be on the outer edge 
of the outermost bright HH-regions in the disk (Figs. [3 
andEJ. Although this definition lacks a specific threshold 
value, the steepness of the radial brightness profile ensures 
that the values for rgF are relatively robust. The resulting 
rsF values are listed in Table 01 

Note that once the correct value of tsf is determined, 
one can use integral (IRAS) far-infrared fluxes or radio 
continuum data-neithcr of which are affected by optical 
thickness-to determine average energy input rates by mas- 
sive SF per unit area, by assuming that all emission arises 
from the star- forming parts of the disks. 

The values for the radii of the radio continuum- 
emitting disks, tcr, were measured at the 5-er confidence 
level of the images by producing cuts through the radio 
continuum emission distribution in the galaxy disks along 
their major axes, as described by us in DLG95. Again, two 
measurements are taken and then averaged. The place- 
ment of the cut through the 1.49 GHz disk emission of 
NGC 4666 is exactly the same as that for the Ha cut 
(Fig. 01. The resulting surface brightness profile as a func- 
tion of galactocentric distance is shown in Fig.|3]as a solid 
line. Details are listed in Table and the resulting mea- 
surements of tcr can be found in Table 
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Table 3. Radii of CR Halos and Comparison with Thin Nonthermal Radio Disks and Ha Disk Emission 
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(1) 


(2) 




(3) 




(4) 




(5) 


(6) 




(7) 


(8) 


Object 


D 




tcr 




fhalo 




rsF 


Ar 










[Mpc] 


["] 


[kpc] 


["] 


[kpc] 


["] 


[kpc] 


[kpc] 


["] 


[kpc] 






r i 

U.KJ 


97& 


1 9 & 




Q 9 


90^ 


Q 3 
y . o 


o . o 


49zL 




99 


NGC 4631 W 


10.0 


492 


23.8 


421 


20.4 


362 


17.6 


6.2 


444 


21.5 


0.67 


NGC 4631 E 




433 


21.0 


338 


16.4 


348 


16.9 


4.1 


444 


21.5 


0.62 


NGC 1808 


10.9 


9±2 


0.47±0.10 


15±4 


0.79±0.21 


8±2 


0.42±0.10 


0.05±0.15 


195 


10.3 


0.0017 


M82 


3.2 


46±3 


0.71±0.03 


60±3 


0.93±0.05 


21±2 


0.33±0.04 


0.38E0.06 


336 


5.2 


0.0039 


NGC 4666 


26.4 


119±4 


15.2±0.5 


90±6 


11.5±0.8 


106±4 


13.6±0.5 


1.6±0.7 


138 


17.7 


0.59 


NGC 4700 


25.5 


63±5 


7.7±0.5 


42±7 


5.2±0.6 


51±7 


6.3±0.5 


1.1±0.6 


90 


11.1 


0.35 


NGC 7090 


11.7 


147±8 


8.3±0.5 


93±12 


5.3±0.7 


74±8 


4.2±0.5 


0.9±0.6 


222 


12.6 


0.35 



The entries in Table |S| are: 

(1) Galaxy name; (2) Distances are based on Ho = 75 km s _1 Mpc -1 and a virgocentric infall velocity of 300 km s" 1 (cf. 
Dahlem et al. 2001); (3) Radial extent of the radio continuum-emitting disk; (4) Radial extent of the radio synchrotron halo; 
(5) Radial extent of the star-forming disk, measured from Ha; (6) Difference between rcn and rsF', (7) Radius of stellar disk 
at 25th magnitude (in the B-band); (8) Scaling factor between star-forming part of the disk, Asf, vs. stellar disk extent, ^25- 



3.2. Radial extent of the radio halos, rhaio 

We measure the radial extent of the radio halos, rhaio, by 
producing cuts through the radio continuum emission dis- 
tribution in the halos, parallel to the galaxies' major axes, 
but offset from the disk planes (DLG95). Measurements 
are again taken at the 5-er confidence level. We show the 
positions of the cuts through the radio halo of NGC 4666 
in Fig. 0] where we follow exactly the same procedure 
as in DLG95 (cf. Fig. 3 in that paper). The placement 
of the cuts is quantified in Table El From the two cuts 



through the halo of each galaxy, one can ideally obtain 
four measurements of rhaio- Because axial symmetry is a 
good assumption for the relatively undisturbed systems 
studied here, the four values can be averaged into one per 
galaxy. Fig. [S] displays, as an example, the points at which 
rsF and rhaio were measured in NGC 4666. In practice, 
some values cannot be used, e.g. because of the presence 
of nearby background sources. Column 6 in Table [3 pro- 
vides the number of measurements used to obtain rhaio for 
each galaxy. The average values for rhaio obtained by us 
(in angular and spatial units) are presented in column 4 




Fig. 3. 12" wide cut through our Ha image (dotted 
line) and through the 1.49 GHz emission distribution of 
NGC 4666 (solid line) along the disk plane (cf. Fig. El, 
demonstrating the measurements of r§F and rcR. Arrows 
mark the positions at which radii were measured; upwards 
arrows at the rcR measurements (at the intersection of the 
solid line with the dotted 5-cr threshold value), downwards 
arrows at the rgp points (at the intersection of the dotted 
data graph with the dashed horizontal threshold line). 



33' 




12 h 45 m 20 s 45 m 10 s 45 m 00 B 

Right Ascension (J2000) 



Fig. 4. Overlay of rectangular boxes with widths of 12" 
for the measurements of rcR and rhaio in NGC 4666 on 
our 1.49 GHz radio continuum image. The central cut is 
the same as in Fig. El the properties of the two parallel 
cuts through the halo are listed in Table El 
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Fig. 5. Cuts through the 1.49 GHz emission distribution 
of NGC 4666 through the radio halo, parallel to the disk 
plane, at z-offsets of ±32" (see Table |5J|. The cut through 
the halo north-west of the plane is shown as a dashed line, 
the south-eastern cut is represented by a solid line. For 
comparison the cut through the Ha image along the disk 
plane is shown again as a dotted line. Arrows mark the 
positions at which radii were measured; upwards arrows 
at the rhaio measurements, downwards arrows at the tsf 
points (same as in Fig.[3J). 



of Tab. |3| Since it is important that the cuts in the halo 
not be contaminated by disk emission (by either beam- 
smearing, a warped disk or a deviation of inclination from 
edge-on [i=90°]), z-offsets where chosen to be more than 
one beam- width. At the same time, the z-offsets must not 
be too large because the signal-to-noise ratio of the data 
drops rapidly (exponentially) away from the galaxy disks. 
Accordingly, cuts were produced at intermediate z-offsets, 
well above the projected disks of the galaxies, but as low as 
possible in order to achieve the highest possible signal-to- 
noise ratios. We list the z-offsets of the cuts from the major 
axes in units of arcsec, kpc and beamwidths in Table |2 

All average values for the three radii (rsF, ^cr an d 
Hiaio) are collated in Tabled The values in the top part of 
the table, separated by a horizontal line, are reproduced 
from DLG95, the lower part of the table contains new 
measurements. 

Our estimates of the uncertainties for the various radii 
are dominated by uncertainties arising from the limited 
angular resolution of the radio images. In addition, the 
observed radio halos are not exactly cylindrical and there- 
fore there is a weak dependence of the radial extent, rhaio ! 
on the z-offset at which the measurement is obtained. 
Because of the constraints mentioned before, this uncer- 
tainty can at present not be avoided. 

Due to signal-to-noise limitations, studies of variations 
of the halo width as a function of z-distance can currently 
not be conducted yet. 



4. Discussion 

Our new results enable us to perform the same studies 
on a sample of five new galaxies as done previously for 
NGC 891 and NGC 4631 (DLG95), with which they can 
be compared: 

1. Check whether there is a relation between the maxi- 
mum radius at which massive SF is observed, tsf, and 
the radial extent of the radio halo, rhaio fSec. 14.1(1 : 

2. calculate the average energy input rate from high- mass 
SF in the disk, , and determine the threshold value 
for the onset of halo emission fSec. 14. 3|) . 

In addition, the slowly increasing size of our sample 
now also allows us to investigate whether the existence 
and properties of gaseous halos depend on a number of 
other global physical galaxy properties fSec. l4.4|l . For this 
study values obtained for other galaxies (DLG95) have 
been reused. 

4.1. Star-forming disk vs. radial halo extent 

As before (DLG95), we find a strong dependence of the 
halo properties on the level of disk activity in our sample 
galaxies. The numbers quantifying these relationships are 
collated in Table|21 A tight and linear correlation between 
the radial extent of radio halos, rhaio, and the areas of ac- 
tive SF in the underlying disks, tsf, is confirmed. This is 
visualised in the left panel of Fig. A weaker and strongly 
nonlinear dependence is visible when plotting rhaio against 
the radial extent of the optical disk at the 25th magnitude 
surface brightness, r25 (right panel of Fig. El . The use of 

is not meaningful in this context, as can be judged by 
considering a study of galaxies with proto-typical circum- 
nuclear starbursts. In these galaxies, the vast difference 
between r^ and tsf (see entries for NGC 1808 and M 82 in 
Table |2J) illustrates clearly that is not a suitable tracer 
of currently ongoing SF activity. This, in turn, supports 
the notion that the observed gaseous halos are created by 
processes related to the young stellar populations (I), not 
the old ones (II). 

The extent of the radio halos, rhaio, of galaxies with 
widespread SF is similar to that of the underlying star- 
forming disks, rgF, because in the outer parts of the galaxy 
disks the SF rates (SFRs) and thereby energy input rates 
are too low to create outflows (cf. DLG95). The rhaio val- 
ues are in most cases slightly smaller than rgF, which is 
probably a sensitivity bias. These data points lie just be- 
low the line of unity in the left panel of Fig. [S] 

A direct comparison of the values of r$F and rhaio 
(Table EJ) suggests that for circumnuclear starbursts 
(NGC 1808, M82) r ha i is slightly larger than r SF . This 
is caused by the fact that the starburst-driven outflow 
cones widen beyond the disk planes in the arche-typical 
hour-glass shape (e.g. Veilleux et al. 2005 and references 
therein). For this reason, classical starbursts will always 
appear slightly above the line of unity in Fig. (left panel). 
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Fig. 6. Left: Radial extent of the star-forming part of the galaxy disks, rgp, against the radial extent of the associated 
halos, rhaio! see Sec. 14. II for details. Right: Radial extent of the optical disk, T25, against fhaio- 



The significance of a direct relation of the halo extent 
and the maximum radius at which a galaxy is actively 
forming stars is that-based on the data used here-halos 
appear to exist only when and where the current massive 
SF rate is high. It corroborates calculations predicting out- 
flows from active SF regions predominantly perpendicular 
to the disk plane due to the pressure gradient along which 
material can flow. At the same time expansion within the 
disk is inhibited by the resistance of swept-up disk mate- 
rial. 

One example of a model making such predictions, al- 
though too simplistic to explain all details, is the "chim- 
ney" model by Norman & Ikeuchi (1989). Any modern 
model of the interstellar medium in galaxies contains a 
halo component, which is created naturally as a conse- 
quence of the buoyancy of disk gas heated by massive 
SF (see, e.g,. Mac Low & Ferrara 1999). Bregman (1980) 
pointed out that the conservation of angular momentum 
during the transport of matter into a galaxy's halo (in 
the form of a "galactic fountain") leads to radial motions, 
thus facilitating an efficient radial re-distribution of metal- 
enriched material within galaxies. 

The major advance of our current observational results 
over DLG95 is that our findings are now firmer, because 
based on significantly more sensitive observations of more 
objects (seven compared to two) than before. 



the radio continuum-emitting part of a galaxy can be ob- 
served: 

1. Diffusion of CRs from sites of massive SF. In the case 
of radial diffusion within the disk, for galaxies in which 
the level of SF drops rapidly at a certain galactocentric 
radius, Ar = r^R — ?"sf is a measure of the radial CR 
diffusion coefficient within the disk, D (DLG95). 

2. The presence of low-level SF in the outermost parts of 
the disks can lead to low surface brightness radio emission. 

Galaxies with widespread SF (in the present sam- 
ple NGC4666, NGC4700 and NGC7090; NGC891 and 
NGC4631 from DLG95) exhibit a significant difference 
between r$F and tqr (see Table |3J). This suggests either 
radial electron diffusion in the disks or low-level SF in the 
outer parts of their disks, or a combination of both. 

On the other hand, the measured differences between 
rgp and tqr within the disks of circumnuclear starbursts 
are insignificant. A likely cause of this similarity could 
be that radial propagation of CRs within the disks is in- 
hibited by the high pressure of the ambient medium and 
possibly by poloidal magnetic fields, which are found in 
powerful starburst galaxies (e.g., Lesch et al. 1990). At 
the same time, there is little or no SF in the disks of 
some circumnuclear starburst galaxies beyond their Inner 
Lindblad Resonance (or turnover radius of rotation); see 
e.g. Combes (1987), Lesch et al. (1990). 



4.2. Extent of the star-forming disk vs. total radio 
continuum-emitting disk 

There are two reasons why a significant difference between 
the extent of an actively star-forming disk and that of 



4.3. Energy input rates into the disk ISM 

Based on rgp we calculate the surface area within the sam- 
ple galaxies over which SF occurs, assuming circular sym- 
metry, Asf — 7r 7g F . This we adopt as the surface area over 
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Table 4. Energy Input by Supernovae Derived from Radio Data 



(1) 




f3l 
\°J 


(4) 


C5l 


f6l 
\y) 


(7) 


\j o j ec L 


IVIOI pn. 


u 


TSF 


-filth 




7Tlt0t 




Type HC3 


[Mpcl 


fkpcl 


[10 21 WHz" 1 ] 


IJ A J 


[10 -3 ergs s _1 cm -2 ] 


NGC 891 


SA(s)b? sp 


9.5 


9.3 


8.2 


0.13 


1.60 


NGC 3044 


SB(s)c? sp 


20.6 


10.0 


6.1 


0.10 


1.02 


NGC 4631 


SB(s)d 


10.0 


17.3 a 


14.7 


0.23 


0.82 


NGC 5907 


SA(s)c: sp 


14.9 


18.1 


2.6 


0.04 


0.13 


NGC 4565 


Sbc 


9.7 


16.9 


1.5 


0.02 


0.09 


NGC 4244 


SA(s)cd? sp 


3.1 


4.5 


0.1 


0.002 


< 0.08 



M82 


I0;Sbrst HII 


3.2 


0.33 


10.1 


0.16 


1550 


NGC 1808 


(R'_l)SAB(s:)b 


10.9 


0.42 


5.9 


0.09 


565 


NGC 4666 


SABc: 


26.4 


13.7 


34.1 


0.53 


3.1 


NGC 3175 


SAB(s)b 


15.9 


3.5 


2.1+ 


0.03 


2.98 


NGC 4700 


SB(s)c? sp; HII 


25.5 


6.3 


2.3 


0.04 


0.99 


NGC 1406 


SB(s)bc: sp 


14.9 


7.9 


3.4+ 


0.05 


0.90 


NGC 3437 


SAB(rs)c: 


25.5 


11.8 


5.7+ 


0.09 


0.69 


NGC 7090 


SBc? sp 


11.7 


4.2 


0.6* 


0.01 


0.59 


NGC 2748 


SAbc HII 


28.7 


14.3 


6.6+ 


0.10 


0.54 


NGC 1421 


SAB(rs)bc: 


31.1 


21.8 


12.8+ 


0.20 


0.45 


NGC 7462 


SB(s)bc? sp 


15.1 


5.4 


0.7+ 


0.01 


0.41 


NGC 1055 


SBb: sp LINER2 


16.0 


18.3 


7.2+ 


0.11 


0.36 


NGC 3717 


SAb: sp HII 


27.1 


23.2 


10.8+ 


0.17 


0.34 


NGC 5297 


SAB(s)c: sp 


34.5 


12.7 


3.3 


0.05 


0.34 


NGC 2613 


SA(s)b 


25.9 


21.2 


5.5 


0.08 


0.21 


NGC 4517 


SA(s)cd: sp 


19.5 


18.4 


1.7 


0.03 


0.08 



The entries in Table 0] are: 

(1) Galaxy name; (2) Morphological type from NED, based on the RC3 catalogue; (3) Distances are based on Ho = 75 km 
s _1 Mpc -1 and a virgocentric infall velocity of 300 km s _1 (cf. Dahlem et al. 2001); (4) Radial extent of the star-forming disk, 
measured from Ha imagery; (5) Total nonthermal radio power at 1.49 GHz (1.384 GHz for galaxies marked with an asterisk 
[*]; 1.43 GHz for galaxies marked with a cross [ + ]); (6) Supernova rate derived from the total radio power, P n th, applying eq. 4 
of DLG95; (7) The energy input per SF area derived from eq. 8 of DLG95, assuming for all galaxies a magnetic field strength 
of 5 /xG. 
Note to Table E| 

a) Average value for both sides of the disk 



which energy is released homogeneously into the ambient 
medium by SF-related process, i.e. stellar winds, super- 
novae and their remnants (Leitherer and Heckman 1995), 
Sgpj. The energy input is assumed to be proportional to 
the nonthermal radio emission and its calculation and nor- 
malisation to unit surface area, E^ = E<$/ A$f, are per- 
formed exactly as described in DLG95. The relevant num- 
bers are collated in Tables |3] and where the top parts 
contain values from DLG95, while the lower parts add new 
measurements to the database. 

Especially for the classical circumnuclear starbursts 
the use of rgF instead of is obviously essential. 
Neglecting this correction would lead to errors in estimates 
of the mean energy input rates per unit surface area, , 
by factors of 1000, typically, as can be seen by inverting 
the numbers in column 8 of Table |3] But also for galaxies 
with widespread SF corrections in E^ by factors of 1.5-5 
occur when using rsF instead of r25 (Tables 01 and . 

Again, there is a clear difference between proto-typical 
starburst galaxies (NGC 1808, M 82) and others with more 
widespread SF. The energy injection rates in circumnu- 



clear starbursts, E^, are orders of magnitude higher than 
in galaxies with widespread SF, leading to much more en- 
ergetic outflows which can be detected at high distances 
above the disk planes, of up to 5-10 kpc (Dahlem et 
al. 1998, Lehnert et al. 1999, Strickland et al. 2004a,b, 
Tiillmann et al. 2006). In these outflows gas can reach 
of order escape velocity (e.g., HAM90) and shock heating 
plays an important role (e.g., Chevalier and Clegg 1985, 
HAM90, Dahlem et al. 1997). 

Outflows from galaxies with widespread SF reach lower 
maximum z-heights of up to 1-5 kpc (Rand 1996, Rossa & 
Dettmar 2000, 2003a,b, Miller & Veilleux 2003), because 
the average energy input rates per unit surface area are 
lower than in starbursts, although the total energy input 
may be similar (see column 5 in Table 0J. 

4.4. Halo dependence on various physical properties 

Several quantities were investigated for potential relation- 
ships. In particular, possible dependences of the halo prop- 
erties on the level of SF in the galaxy disks, as measured 
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Table 5. Fundamental Properties of Galaxies Shown in Figures GUIS 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Object 


/60 


/ioo 


/ioo 


Lfir 




rtiK 


M K 


M K /A 25 




[Jy] 


[Jy] 


[10 10 Lq] 


[10 10 L ] 


[mag] 


[10 10 Mq] 


[Mq pc" 2 ] 



NGC 891 


61.10 


198.63 


0.31 


1.27 


1.44 


5.94 


3.2 


27.0 


NGC 1055 


23.27 


60.09 


0.39 


1.21 


1.18 


7.15 


3.0 


30.6 


NGC 1406 


11.91 


27.00 


0.44 


0.51 


0.54 


8.61 


0.7 


31.9 


NGC 1421 


11.20 


24.30 


0.46 


2.03 


3.43 


8.40 


3.6 


44.5 


NGC 1808 


87.81 a 


137.20 a 


0.64 a 


1.70 a 


0.70 


6.66 


2.2 


65.9 


NGC 2613 


6.31 


25.42 


0.25 


1.10 


6.37 


6.83 


10.6 


47.6 


NGC 2748 


7.95 


19.44 


0.41 


1.30 


1.70 


10.32 


0.5 


10.6 


M82 


1271.32 


1351.09 


0.94 


1.87 


0.26 


4.67 


1.2 


138.9 


NGC 3044 


10.47 


21.16 


0.49 


0.81 


1.34 


8.98 


0.9 


15.6 


NGC 3175 


13.07 


28.20 


0.46 


0.62 


0.73 


7.79 


1.6 


61.6 


NGC 3437 


12.15 


20.62 


0.59 


1.33 


0.92 


8.88 


1.5 


57.6 


NGC 3717 


10.52 


23.84 


0.44 


1.48 


2.67 


7.52 


6.1 


35.0 


NGC 4244 


4.20 


16.06 


0.26 


0.01 


0.17 


7.72 


0.07 


4.8 


NGC 4517 


6.92 


20.20 


0.34 


0.57 


3.61 


7.33 


3.8 


14.3 


NGC 4565 


9.83 


47.23 


0.21 


0.27 


1.92 


6.06 


3.0 


22.0 


NGC 4631 


82.90 


208.66 


0.40 


1.67 


3.21 


6.67 


2.2 


15.2 


NGC 4666 


37.34 


82.88 


0.45 


4.92 


3.32 


7.06 


8.9 


91.1 


NGC 4700 


3.05 


5.36 


0.57 


0.34 


1.54 


9.78 


0.7 


17.4 


NGC 5297 


1.81 


7.99 


0.23 


0.59 


2.95 


9.89 


1.1 


5.4 


NGC 5907 


12.02 


42.37 


0.28 


0.64 


2.61 


6.76 


3.7 


18.4 


NGC 7090 


5.88 


17.97 


0.33 


0.18 


0.97 


8.16 


0.6 


12.7 


NGC 7462 


3.11 


6.32 


0.49 


0.13 


0.55 


9.65 


0.3 


10.0 



The entries in Tableware: 

(1) Galaxy name; (2) IRAS 60 ^m band flux density; from NED; (3) IRAS 100 /im band flux density; from NED; (4) /eo/Zioo 
flux ratio; (5) Far- infrared luminosity, based on the standard formulas used by us earlier (Dahlem et al. 2001); (6) Optical 
blue luminosity, determined as Lb = 89.13 x io-°- 4xM b i wnere t he absolute blue magnitude is the apparent blue magnitude, 
K-corrected, corrected for Galactic extinction and scaled by the distance of the object, Mb = m,g -25-5xlog(Z?); values 
are from NED. This definition provides an estimate of the blue luminosity (isL„) at 4400 A; (7) Apparent total magnitude in 
the K-band (centered at 2.17 /im) from 2MASS; (8) Mass in galactic disk derived from the K-band magnitude, as described in 
Sec. I4.4.T1 (9) Average mass surface density within the optical (r2s) disk. 
Note to Tabled 

a) Not used for studies of FIR properties due to unknown contribution from a low-luminosity AGN. 



by the energy deposition rate E^, and on galaxy mass 
surface density were tested. 

4.4.1. Energy injection rate vs. galaxy mass surface 
density 

Lower-mass galaxies not only have lower escape velocities 
than more massive ones, but the onset of disk-halo in- 
teractions will also be easier for a shallower gravitational 
potential. Therefore, the existence of a radio halo is ex- 
pected to depend, apart from the energy input from the 
disk (see above), on the gravitational force that particles 
must overcome to escape from a galaxy disk. The grav- 
itational potential that particles experience when being 
lifted above the disk is proportional to the mass surface 
density in the disk. We approximate the mass in the disk 
by the stellar mass derived from near-infrared emission, 
which traces old stars that are responsible for the bulk of 
the stellar mass. 



The stellar mass is deduced from the total K-band 
magnitude, mx, from the Two Micron All Sky Survey 
(2MASS). The luminosity in this band, Lk, was calcu- 
lated as 

L K = 2.8 x 10 10 p^lO- 4mK [L ], (1) 

which gives an estimate of the luminosity (y L u ) at 2.17 
/mi. 

From the luminosity the stellar mass was derived as- 
suming optically thin emission (which may lead to a slight 
underestimate of the NIR emissivity) and a mass-to-light 
ratio of 3, which is indicative of a stellar population of 
age 5.9 Gyr (Bruzual & Chariot 2003). The K-band mag- 
nitudes and Mk values for the galaxies in our sample are 
listed, together with other parameters, in Table 

The significance of the rate at which energy is de- 
posited into the ISM of spiral galaxies, E^, on the ini- 
tiation of disk-halo outflows, and thereby the creation of 
gaseous halos, has been described above (Sec. I4.3JI . Only 
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Fig. 7. Left: Energy input per unit star- forming area, Asf, as a function of stellar mass surface density, as derived 
from K-band imagery. Filled symbols indicate galaxies with radio halos, open symbols are galaxies without radio 
halo detections. The dashed line roughly indicates the transition zone between galaxies with and without radio halos. 
Right: The same as on the left, but using A25 instead of Asf as an approximation of the area over which energy is 
deposited into the disk ISM. 



galaxies with energy input rates above a certain thresh- 
old value can start disk- halo interactions. Note that E^ 
values can be converted directly into the commonly used, 
but more indirect, SF rates per unit area, which are often 
measured in units of M yr _1 kpc~ 2 . 

In order to investigate the influence of both the galaxy 
mass surface density and the energy input rates, we plot in 
the left panel of Fig.[7|the stellar surface density, M-^/ A25, 
where A25 = 7rr| 5 is the area of the optical disk, versus 
Ej^ 1 for the galaxies in our extended sample. This plot is 
similar to the one shown in Lisenfeld et al. (2004), however 
now using the stellar mass in the disk, rather than the 
total dynamical mass. In addition, the present plot also 
includes more new galaxies from the present paper. 

As expected, we note a clear division between galax- 
ies with and without radio halo: Galaxies with low mass 
surface densities and high energy input rates have radio 
halos, as opposed to galaxies with high mass surface den- 
sities and low energy input rates. The line drawn indicates 
the transition between both regimes. 

The significance of the observed dichotomy is that, be- 
cause of the smaller force that particles must overcome and 
the associated lower escape velocities, it is easier for low- 
mass galaxies to have gaseous halos around them than 
massive ones (e.g. Heckman et al. 1995, Martin 1997, 
1998). Once the energy input rate is high enough, also 
spiral galaxies with a high mass surface density can cre- 
ate and maintain gaseous halos. The energy input levels in 
classical starburst galaxies are shown to be so high (Fig.0 



left panel) that gravitational forces are negligible, as al- 
ready assumed by Chevalier and Clegg (1985). 

These new findings refine our statement (DLG95) that 
there is "a threshold energy input rate" above which out- 
flows can occur by including the galaxy mass surface den- 
sity as an additional parameter. The energy threshold 
value determined here for galaxies of given mass surface 
density for the onset of disk-halo interactions can be used 
in the future to predict whether or not a certain galaxy 
can be expected to have a gaseous halo, i.e. whether the 
"breakout condition" , as defined by Norman & Ikeuchi 
(1989), required for outflows to occur is fulfilled. The de- 
tected E^ t vs. Mk/^25 relation is reliable, because it does 
not depend strongly on the adopted galaxy distances, both 
quantities essentially being derived from surface bright- 
nesses. 

In this context it is important to note that using r25 
instead of tsf to measure E^ would lead to the displace- 
ment of objects along the x-axis to lower values of E^, 
which would skew the detected E^ vs. AIkA^s) rela- 
tionship. This is illustrated in the right panel of Fig. [7| 
Accordingly, no such clear separation between galaxies 
with and without a halo is visible there and the place- 
ment of data points for classical starbursts is completely 
wrong. 

Strongly interacting galaxies must be excluded from 
studies of the dependence of halo properties on the un- 
derlying SF activity, because of the potential influence of 
the gravitational force of a companion galaxy on the gas 
kinematics and thus disk- halo interactions. At the present 
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time we do not have sufficient information to disentangle 
the two effects. 

4.4.2. Localised outflows vs. large-scale halos 

The values of E^ derived above are averages. These are 
a measure of the global energy input rates into a galaxy's 
ISM required to form large-scale halos. 

As calculated by us earlier (DLG95), local values of 
E^ lt in active SF regions are more likely to be of order 
10~ 2 ergs s _1 cm~ 2 . Fig. [7| suggests that local outflows 
above giant H II regions, with E^ ~ 10~ 2 ergs s _1 cm~ 2 , 
should always be possible, even in otherwise quiescent 
galaxies with high mass surface densities. This explains 
why a galaxy such as NGC 4565, with a low SF rate and 
thereby low average energy input level (see Fig. UJ, can 
exhibit a small number of dust filaments perpendicular to 
its disk plane which might be tracers of localised outflows 
of material. Note that, apart from the non-detection of a 
radio halo, NGC 4565 is also deficient in Ha and [OIII] in 
the disk- halo interface (Rand et al. 1992 and Robitaille et 
al. 2006), which again argues against the existence of a 
large-scale gaseous halo. 

The decisive difference between a circumnuclear star- 
burst and an individual localised outflow in a disk with 
more widespread SF is the total amount of energy injected 
into the ISM. The energy injection rates, E^ , are yet 
higher in starbursts than in Giant Extragalactic H II re- 
gions (GEHRs; Kennicutt 1984). Also, since classical star- 
bursts comprise a number of star-forming regions each of 
which equals a GEHR, energy (and mass) is deposited 
into the ambient medium over a longer time than in an 



individual GEHR. Hence, while for an individual GEHR 
the total energy can amount to of order 10 54 ergs, that of 
a classical starburst is orders of magnitude larger, nearer 
10 56 ergs. In the case of ultra-luminous infrared galaxies 
(ULIRGs) up to 10 58 " 59 ergs are deposited (HAM90). 

4.4.3. Compactness of the star formation distribution 

The left panel of Fig. [5J suggests that there is no clear link 

between the far-infrared luminosities, Lfir, and the ratio 

of the surface areas of the star-forming parts vs. the total 

extent to 25th optical blue magnitude of the galaxy disks, 

Asf/^25, in our sample galaxies. In fact, when comparing 

the position of M82 on one hand with those of NGC 891 

and NGC 4631 on the other, one can see that for widely 

different SF distributions galaxies can have the same Lfir- 
i 

However, in the right panel of Fig. [SJ although again 
there is no clear correlation between the two quantities, 
a trend is indicated that galaxies with radio halos have 
smaller A^p/A^^ values and at the same time higher 
Lfir/Lb ratios than the ones without halos. This sug- 
gests that a compact SF distribution favours the creation 
of a gaseous halo. 



1 NGC 1808 was excluded from this part of the investigation, 
because the contribution of its low-luminosity Seyfert 2 AGN 
to the total FIR luminosity is unknown. In the determination of 
its radio flux density (above) the flux contribution of the AGN 
was subtracted by fitting the central peak using the AIPS task 
JMFIT. 
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Fig. 9. Left: Agp/t^bs plotted vs. the IRAS 60 (im to 100 /im flux ratio, which is a measure of the average warm 
dust temperature. Right: Energy input rates per unit surface area, Ejg*, plotted vs. feo/fioo- NGC 1808 is not shown, 
because the contributions of its AGN and starburst to the FIR flux densities could not be separated. 



The left panel of Fig. [5] shows another trend, namely 
for galaxies with radio halos to have both compact SF 
distributions and warm average dust temperatures (as re- 
flected by their high IRAS 60 /im to 100 /im flux ratios, 
/eo/Aoo)- The right panel of Fig. [§| which is similar to 
Figs. 38 and 39 by Dahlem et al. (2001) and Fig. 5 by 
Rossa and Dettmar (2003a), indicates that galaxies with 
radio halos also have high energy input rates into their 
ISM and correspondingly high average dust temperatures. 

These relationships show that any combination of high 
-£'FIr/-E'B ratio, compact SF distribution and high feo/ fwo 
flux ratio ( ~ 0.4) is a good tracer of galaxies with radio 
halos. This explains why high FIR flux densities and high 
feo/ i too flux ratios (indicating the presence of warm dust) 
are such good selection criteria to find, after rejection of 
AGNs, galaxies with SF-driven gaseous halos (Dahlem et 
al. 2001). 

4.4.4. Energy input rate vs. halo z-height 

A search for a dependence of the radio halo z-scale heights 
on the level of energy input in the disk did not yield con- 
clusive results. Since the uncertainties involved are yet too 
large to make firm statements, this topic will not pursued 
any further at this time. 

4.5. Future work 

It is unlikely that large numbers of new targets fulfill- 
ing the stringent selection criteria for studies like this 
will be found in the near future. The main search criteria 



(FIR flux and colour) become more unreliable when go- 
ing to more distant and thus fainter sources. In addition, 
the detection and proper imaging of gaseous halos in dis- 
tant galaxies is extremely time-consuming on the present 
generation of telescopes. A much larger sample will be- 
come available only through the advent of the next gener- 
ation of more powerful telescopes (both optical and radio). 
Progress in this field at the present time can be made by 
obtaining more complete databases for the galaxies pre- 
sented here (and a few others known already, including e.g. 
NGC 253) and by finding ways to measure halo properties 
in face-on (or at least low inclination) galaxies. 
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